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Abstract 

It was recently shown (Keating & Prado, Proc. R. Soc. Lond. A 457, 1855-1872 
(2001)) that, in the semiclassical limit, the scarring of quantum eigenfunctions by 
classical periodic orbits in chaotic systems may be dramatically enhanced when the 
orbits in question undergo bifurcation. Specifically, a bifurcating orbit gives rise to 
a scar with an amplitude that scales as h a and a width that scales as ft , where a 
and u are bifurcation-dependent scar exponents whose values are typically smaller 
than those (a = u = 1/2) associated with isolated and unstable periodic orbits. We 
here analyze the influence of bifurcations on the autocorrelation function of quantum 
eigenstates, averaged with respect to energy. It is shown that the length-scale of the 
correlations around a bifurcating orbit scales semiclassically as h l ~ a , where a is the 
corresponding scar amplitude exponent. This imprint of bifurcations on quantum 
autocorrelations is illustrated by numerical computations for a family of perturbed 
cat maps. 



* Short title: Bifurcations and autocorrelations 
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1 Introduction 



The local statistical properties of quantum eigenfunctions VVi(q) in classically chaotic 
systems are often modelled in the semiclassical limit by treating the eigenfunctions as 
random superpositions of plane waves (Berry 1977). This allows, for example, for value 
distributions and autocorrelation functions to be computed. However, it ignores contri- 
butions from the scarring of the eigenfunctions by classical periodic orbits (Heller 1984, 
Kaplan 1999). Understanding these contributions has been one of the central goals in 
quantum chaology. Our purpose here is to study the influence of bifurcating periodic 
orbits on the autocorrelation of eigenfunctions. 

The semiclassical theory describing the scarring of wavefunctions by isolated and un- 
stable periodic orbits was developed by Bogomolny (1988), who derived an explicit formula 
for the contribution from a given such orbit to |^ n (q)| 2 , averaged locally with respect to 
position and energy. Berry (1989) extended this approach to the Wigner functions associ- 
ated with energy eigenstates, again averaged locally with respect to the energy (position 
averaging is unnecessary in this case). Berry's formula also describes the influence of 
isolated and unstable periodic orbits on the autocorrelation of the eigenfunctions, be- 
cause the Wigner function and the autocorrelation function {ip n (q, + L)^*(q)) are fourier 
transforms of each other (explicit formulae my be found in Li & Rouben 2001). It was 
shown in Keating (1991) that Bogomolny 's and Berry's formulae are exact, rather than 
semiclassical approximations, for quantum cat maps. 

The extension of Bogomolny's theory to bifurcating periodic orbits was made in Keat- 
ing & Prado (2001). There it was shown that orbits undergoing bifurcation give rise 
to scars in |-0 n (q)| 2 (averaged locally with respect to position and energy) with a width 
that scales as h w and an amplitude that scales as h a when h — > 0, where u and a are 
bifurcation dependent scar exponents whose values are positive and less than or equal 
to those (cu — a — 1/2) corresponding to isolated and unstable orbits. For most of the 
generic bifurcations u < 1/2 and a < 1/2, and so bifurcations may be said to give rise 
to superscars. This mirrors the dominant influence they have previously been shown to 
exert on spectral fluctuation statistics (Berry et al. 1998, Berry et al. 2000). 

The question we address here is how bifurcations influence the autocorrelation of the 
eigenfunctions. Our main result is that the length-scale of the correlations around a 
bifurcating orbit scale as ti 1 ~ a when h—>0, where a is the scar amplitude exponent. By 
comparison, the random wave model predicts correlations with a length-scale of the order 
of the de Broglie wavelength (i.e. of the order of h). We illustrate the influence of scarring 
on the quantum autocorrelation function with numerical computations for a family of 
nonlinearly perturbed cat maps (Basilio de Matos & Ozorio de Almeida 1995; Boasman 
& Keating 1995). 

The bifurcation of periodic orbits is a characteristic phenomenon in systems whose 
phase space is mixed. It has been conjectured by Percival (1973) that the energy eigen- 
states of mixed systems separate semiclassically into those associated with stable islands 
(regular) and those associated with the chaotic sea (irregular). The random wave model 
has recently been extended to irregular states in mixed systems by Backer & Schubert 
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(2002) (BS1 in the references). The autocorrelation of these states was investigated in this 
context in Backer & Schubert (2002) (BS2 in the references). Our results, concerning the 
length scale of the correlations, are expected to apply to the autocorrelation of irregular 
eigenstates in the vicinity of superscars. 



2 Autocorrelation Formulae 

Our purpose here is to derive semiclassical autocorrelation formulae for bifurcating peri- 
odic orbits which generalize those obtained by Berry (Berry 1989) (see also Li & Rouben 
2001) for periodic orbits far from bifurcation. 

The autocorrelation function for the eigenfunction ^ n (q) of a given quantum Hamil- 
tonian corresponding to an energy level E n is 

c n (L) = ^ n (q+L)^(q)>, (1) 

where (...) denotes an average over position q. We shall be particularly concerned with 
the case when the autocorrelation function is averaged over a group of eigenstates with 
energy levels near to E, that is with 

C(L, E,e) = C "( L K ( E ~ E n) , (2) 

n 

where 8 £ (x) is a normalized, Lorentzian-smoothed 5-function of width e. The right-hand 
side of (j2J thus corresponds to a sum over eigenstates for which E n lies within a range 
of size of the order of e centred on E. Semiclassically, it is approximately the average of 
c n (L) multiplied by d(E), the mean level density, which, for systems with two-degrees-of- 
freedom, is given asymptotically by 

as h — > 0, where 

V(E) = [d(E- H(p, q))dV 2 p (4) 



and H (p, q) is the classical Hamiltonian. 

Our approach follows the one developed by Keating & Prado (2001) for the special 
case L = 0. It is based on the identity 

^.t „ m ^ ^ n (q + L)V£(q) 

G(q + L, q; E) = > , (5) 

where G is the Green function. This implies that 

E + L)C(q)4 (E - E n ) = ~ {G(q + L, q; E + is) 

n 

-{G(q,q + L;E + ie)T}. (6) 



3 



For simplicity, we shall concentrate here on systems with time-reversal symmetry. In this 
case flB} reduces to 



V ^ n (q + L)^(q)* e (£ - S n ) = -- Im G(q + L, q; £ + ie). (7) 
n 

(As will become apparent, our analysis generalizes straightforwardly to non-time-reversal- 
symmetric systems.) The connection with classical mechanics is achieved using the semi- 
classical approximation to the Green function. For systems with two-degrees-of-freedom, 
this is 

G(q + L, q; E) « -== V J\D^\exp \ ^S 7 (q + L, q; E) - ^ I, (8) 

ZftV 27TZAI "~ v L " ^ J 

where the sum includes all classical trajectories from q to q + L at energy E, 5 7 is the 
action along the trajectory labelled 7, 

( d 2 S 1 d 2 S 1 \ 
K£ (9) 
dEdq dE 2 J 

with q' = q + L, and /x 7 is the Maslov index (Gutzwiller 1990). 

The autocorrelation formula we seek follows from averaging (J2J) with respect to q. On 
the left-hand side this gives C(L, E,e). On the right-hand side, the q-average (via the 
stationary phase condition) semiclassically selects from the classical orbits with initial and 
final positions differing by L - henceforth we call such trajectories L-trajectories - those 
for which the initial and final momenta are equal. Our main concern will be with cases 
when |L| vanishes in the semiclassical limit. Then the classical orbits contributing to the 
q-average are close to periodic orbits. To first approximation one can thus find them 
by linearizing about the periodic orbits. Essentially, when a periodic orbit is isolated, 
this corresponds to expanding the action up to terms that are quadratic in the distance 
from it. Specifically, if M p is the monodromy matrix associated with a periodic orbit p, 
then in a section transverse to the orbit, in which the coordinate is y and the conjugate 
momentum is p y (we assume here two degrees of freedom, for simplicity), 

where L y is the component of L in the direction of the y-axis. For a given choice of 
L y , this allows one to determine the coordinates of the associated orbit in the vicinity 
of p which contributes to C(L,E,e). Thus one may view the autocorrelation function 
semiclassically as a sum of contributions from the periodic orbits. The contribution from 
each periodic orbit is, within the linear approximation, stationary in L y when p y = 
(i.e. on the periodic orbit itself), that is when L y = 0. This may be seen explicitly 
by expanding the action in (jJJ) up to quadratic terms around the orbit p. Let z be 
the coordinate along p. Here and in the subsequent analysis we assume that periodic 
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orbit contributions to the autocorrelation function are ^-independent (this is correct in 
determining the h dependence of these contributions when |L| — > as h — > 0, because 
the ^-dependence comes in only via the classical dynamics, except at self-focal points, 
which do not affect the leading order semiclassical asymptotics of the q-average). The 
q-average may then be computed (it is a Gaussian integral) , giving the contribution from 
p to be proportional to exp(ii] p Ly/h), where r\ v is determined by M p . It is clear that 
this is stationary (i.e. centred) at L y = 0. More importantly for us, it also shows that 
the contribution from an isolated periodic orbit has a length-scale that is of the order of 

(in contrast to the random wave model component, for which the length-scale is of 
the order of K). This picture coincides precisely with Berry's (Berry 1989) for Wigner 
functions. The question we now address is what is the length-scale of the contribution to 
the autocorrelation function from bifurcating (i.e. non-isolated) periodic orbits. 

One might imagine that the autocorrelation formulae for bifurcating orbits could be 
obtained easily by expanding the action in ((HI) to higher order than quadratic. This is 
true for some of the simpler bifurcations (e.g. the codimension-one bifurcations of orbits 
with repetition numbers r = 1 and r = 2 - see the examples in Section 3), but for 
more complicated bifurcations it is incorrect, because for these the linearized map 
is equal to the identity, which cannot be generated by the action S^q', q; E) . Thus it 
is difficult to build into the semiclassical expression for the qq'-representation of the 
Green function a well-behaved description of the nonlinear dynamics which the linearized 
map approximates. The solution to this problem, found by Ozorio de Almeida & Hannay 
(1987), is to transform the Green function to a mixed position-momentum representation, 
and this is the approach we now take. 

The method just outlined involves, first, Fourier transforming G(q', q; E) with re- 
spect to q'. This gives the Green function in the q p'-representation, G(p', q; E) (p' is 
the momentum conjugate to q'). The semiclassical approximation to G takes the same 
form as (JHJ), except that S(q', q; E) is replaced by the q p'-generating function S(p', q; E). 
G(q', q; E) may then be rewritten, semiclassically, as the Fourier transform of this expres- 
sion with respect to p'. The result, for the semiclassical contribution to G(q + L, q; E) 
from classical orbits in the neighbourhood of a bifurcating periodic orbit, takes the fol- 
lowing form. 

Consider the general case of a codimension-K bifurcation of a periodic orbit with rep- 
etition number r. As before, let z be a coordinate along the orbit at bifurcation, let y be a 
coordinate transverse to it, and let p y be the momentum conjugate to y, so that y and p y 
are local surface of section coordinates. We again assume the semiclassical ^-scaling of the 
length-scale of the contribution from the bifurcating orbit to the autocorrelation function 
to be ^-independent. Let Q rt K(y,P y , x) be the normal form which corresponds to the local 
(reduced) generating function in the neighbourhood of the bifurcation (Arnold 1978; Ozo- 
rio de Almeida 1988), where x = (x l5 x 2 , ■ ■ ■ , xk) are parameters controlling the unfolding 
of the bifurcation. Then, up to irrelevant factors, the contribution to G(q + L, q; E) is 



(as already stated, we are here interested in determining the ^-dependence of the length- 




(11) 
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scale of contributions to the autocorrelation function, and so have neglected terms in JTTJ), 
such as an ^-independent factor in the integrand, which do not influence this). Note that it 
follows from the definition of the generating function that the stationary phase condition 
with respect to p y selects classical orbits for which the initial and final ^-coordinates 
differ by L y ; that then the stationary phase condition with respect to y (which gives 
the semiclassical contribution to the ^/-average needed to compute the autocorrelation 
formula) further selects from these orbits those whose initial and final momenta, p y , are 
identical; and that the orbits contributing at the point where G r% K is stationary with 
respect to L y are those for which p y = 0. 

We now follow the steps taken in Keating & Prado (2001) for the case when L y = 0: 
first, rescale y and p y to remove the 1/h factor from the dominant term (germ) of $> r ,K in 
the exponent, and then apply a compensating rescaling of the parameters Xi,X2, • • • ,%k 
to remove the ^-dependence from the other terms which do not vanish as h — > 0. Finally, 
we then rescale L y to remove the ^-dependence of the second term in the exponent in 
(|TI|) . Crucially, this term only involves p y , and so if in the initial rescalings 

Py -> H a ^ Py (12) 

then the length scale of the corresponding contribution to the autocorrelation function is 
of the order of h l ~ ar ' K . This is our main general result. It is important to note that a rt x 
is precisely the exponent found for the scarring amplitude in Keating & Prado (2001). 
This then implies a sum rule for the scaling exponents of scar amplitudes and the length- 
scale exponents of wave function correlations reminiscent of those found in the theory of 
critical phenomena. It is worth emphasizing that the correlation length-scale is related 
to the amplitude exponent, rather than the width exponent (denoted u r> K in Keating & 
Prado 2001). 

Values of a r> K were computed by Keating & Prado (2001) for the generic bifurcations 
found in two-degree-of-freedom systems using the appropriate normal forms. For conve- 
nience, we list here the results. First, for the r = 1 (i.e. saddle-node) bifurcations which 
correspond to cuspoid (i.e. corank 1) catastrophes, a^K = \- Values of a for the generic 
bifurcations with r > 1 and K — 1 are listed in Table 1, and those for r > 1 and K = 2 
in Table 2. Finally, a r> K = 2 (k+i) m the g enera l case of bifurcations of orbits for which 
r > 2/\ • 2. 



r 


2 


3 


> 4 




1/2 


1/3 


1/4 



Table 1: Values of av,i for the generic, codimension-1 bifurcations. 

Note that in all cases listed above < a < 1/2, and that in most cases a < 1/2. Thus 
the length scale of the fringes h l ~ ar > K is always semiclassically equal to or smaller than 
for non-bifurcating orbits, and in most cases is smaller. In all cases it is semiclassically 
larger than the de Broglie wavelength. 
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r 


2 


3 


4 


5 


> 6 


«r,2 


1/2 


1/4 


1/4 


1/5 


1/6 



Table 2: Values of a r ,2 for the generic, codimension-2 bifurcations. 
The rescaling of y in (fTTl) was defined in Keating & Prado (2001) by 

y - tf^y, (13) 

where w r) K is the scar width exponent. In the case of the autocorrelation function, 
one must average over y. This leads to the amplitude of G t> k being proportional to 
ffXr,K+Ur,K-i yy e remar k that the exponent of h in the amplitude is related to ampli- 
tude exponent (3 (Berry et al. 2000) of the fluctuations in the spectral counting function 
associated with the bifurcation in question, because 

Oi r ,K + ^r,K - 2 = -{fi r ,K + 1) (14) 

(see equation (2.25) in Keating & Prado 2001). 

The rescaling of x n in (fTll was defined in Keating & Prado (2001) by 

x n -> h a ^ K x n . (15) 

Here, the exponents a describe the range of influence of the bifurcation in the different 
unfolding directions x n , and their sum 

K 

7r,K = ^2 CTn,r,K (16) 
n=l 

describes the /i-scaling of the -fT-dimensional x-space hypervolume affected by the bifur- 
cation. 

3 Perturbed cat maps 

We now illustrate some of the ideas described in the previous section by focusing on a 
particular example: a family of perturbed cat maps. The influence of bifurcations on the 
scarring of eigenfunctions was explored numerically for these systems in Keating & Prado 
(2001). (The influence of tangent bifurcations on localization in wavefunctions was also 
investigated numerically for a different family of maps by Varga et al. 1999). We focus 
here on eigenfunction autocorrelations. 
The maps we consider are of the form 
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where S' p (q) is the first derivative of a periodic function with period 1, k is a parameter 
determining the size of the perturbation, and q and p are coordinates on the unit two- 
torus which are taken to be a position and its conjugate momentum. We will consider 
two cases: S p (q) = sin (2iiq) and S p (q) = cos(2nq). 

These maps are Anosov systems. In both cases, for k < ft max = (v3 — 1)/Vu ~ 0.33 
they are completely hyperbolic and their orbits are conjugate to those of the map with 
k = (i.e. there are no bifurcations). Outside this range, bifurcations occur, stable islands 
are created, and the dynamics becomes mixed (Berry et al. 1998, Keating & Prado 2001). 

The quantization of maps like (fTTj) was developed by Hannay & Berry (1980), when 
k = 0, and Basilio de Matos & Ozorio de Almeida (1995) for non-zero k. The quantum 
kinematics associated with a phase space that has the topology of a two-torus restricts 
Planck's constant to taking inverse integer values. The integer in question, N, is the di- 
mension of the Hilbert space of admissible wavefunctions. With doubly periodic boundary 
conditions (see, for example, Keating et al. 1999), these wavefunctions in their position 
representation have support at points q = Q/N, where Q takes integer values between 
1 and N. They may thus be represented by A^-vectors with complex components. The 
quantum dynamics is then generated by an iV x unitary matrix U whose action on the 
wavefunctions reduces to (fTTj) in the classical limit; for example 



u Q2,Qi = -7^7 exp 



2ni iN 

-TriQl - Q1Q2 + Ql) + sin^Qx/N + tti//2) 



(18) 



where v — gives the sine perturbation and v — 1 the cosine perturbation. This matrix 
plays the role of the Green function of the time-dependent Schrodinger equation for flows. 

Denoting the eigenvalues of U by e lBn and the corresponding eigenfunctions by ty n (Q), 
we have, using the fact that the map (fTTjl is invariant under time-reversal and hence that 
U is symmetric, that 

N N 00 

C(L, 9,e) = J2 Cn(L)5 £ (6 - 6 n ) = 5 L>0 + ^ReJ2Y. U Q+l,q ™P H** - ek) (19) 



N 

n=l Q=l k=l 



where 



1 - 

Cn(L) = -J2^(Q + L)%(Q) (20) 

Q=l 

is the eigenvector autocorrelation function and 

1 - e~ £ cosx 

O e (X) — n 21 

v ' l + e- 2e -2e- £ cosx v ' 

is a periodized, Lorentzian-smoothed 5-function of width e (Keating 1991). Equation 
(jTHjl is the analogue for quantum maps of jZJ). C(L,9,e) corresponds, approximately, to 
times the local n- average (over a range of size of the order of e) of c n (L). For e large 
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enough, the dominant contribution to (TH^I comes from the k — 1 term in the sum on the 
right. We thus define 



1 N 

(L, 9, e) = — U Q+LiQ exp (-id - e) 



(22) 



Q=l 



and, for computational purposes, 

1 N 

c r (L, 9, e) = — Re ^ u Q+l,q exp {-id - e) 



(23) 



We may substitute (fT8| directly into (|22| . In the semiclassical limit, as N — > oo, the 
Q-average selects regions close to stationary points of the phase of (fTHl) . When L = 
these stationary points coincide with the positions of the fixed points of the classical map 
(fTTj) : otherwise the stationary solutions correspond to /-trajectories - in this case classical 
trajectories with initial and final coordinates separated by / = L/N and with initial and 
final momenta equal to each other. These satisfy 



K , 7T, 

7 — t COS (ZTTq- + V— 

2vr v H3 2- 



(24) 



for integers j such that < qj(l) < 1 (see, for example, Boasman & Keating 1995). When 
/ — > (which is the most important regime, given that we wish to focus on correlation 
length-scales that vanish semiclassically) the /-trajectories obviously lie close to fixed 
points of the map. 

Expanding the phase of (jT£l) around qj up to cubic terms gives 



U, 



Q+L,Q 



exp 



27T/.Y,S', -i\ ( 2 - «sin [2^ + u- 

2n 2 iN 



71 



-Kcos^nqj + u-jy 



(2b) 



where 



y = q-q j {l) (26) 

and 2nNSj denotes the phase evaluated at qj(l). Provided that 2 — k sm{2ixqj + i>ix /2) ^ 0, 
this approximation is dominated by the quadratic term in the exponent when y is small. It 
thus describes complex-Gaussian fringes around the classical /-trajectories with a length- 
scale (in terms of y) of the order of N~ l l 2 . 

As discussed in Section 2, the centres of the peaks of the autocorrelation function are 
given by contributing /-trajectories that start at q, end at q + 1, and satisfy p' = p = in 
the first iterate of the map (|T7|l . From (fTT|) . the peak centres are thus solutions of 

l = q^ cos(2vrg + vtt/2) - n (27) 
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and 



3g + - cos (2vrg + z/vr/2) - m = (28) 

7T 

where m and n are integers. For the linear mapping, when k = 0, the peaks are centred 
at I — q — to/3 with m = 0, 1, 2. These peaks may be seen in Figure 1, where the dashed 
line represents an evaluation of (|2~3|) and the bold line a convolution of the data with a 
normalized Gaussian of width 0.007. 

For k > the peak location equations, ([27]) and (|28|) . can be solved numerically. As 
k increases bifurcations of periodic orbits take place. It is important at this stage to 
distinguish between the birth of new periodic orbits as the perturbation parameter k is 
varied (the usual bifurcation) and the birth of new /-trajectories as / varies (a process 
which we call an /-caustic). As / — > 0, the /-trajectories are close to periodic orbits, 
/-caustics coincide with periodic orbit bifurcations in the limit. 

Our main purpose now is to illustrate quantitatively the influence of bifurcating pe- 
riodic orbits on the autocorrelation function. We shall study in detail two bifurcations: 
one a first order (tangent) bifurcation, and the other a second order bifurcation. In order 
to do so, we will treat separately the sine and cosine perturbations. 



3.1 Sine perturbation v = 

For k > 0, the positions of the peaks in the autocorrelation function deviate from those 
(0, 1/3 and 2/3) of the unperturbed map (k = 0). The real part of the autocorrelation 
function for k = 1 is shown in Figure 2 (a). The peak close to / = comes from 
the contribution of the period-1 fixed point j = in (|2~i|) . Using a stationary phase 
approximation with respect to L in (|25|) gives the position of the peak to be 



k cos (27ro (0)) , „ ,„ . 

/ « v , y 1 NN mod 1 29 

tt3 - /tsin(27rg (0)) V ; 

in the limit as / — > 0. For k sin(27rg.,)) < 2, the second order term in the Taylor expansion 
(|2~5|) dominates and the semiclassical formula for the autocorrelation function is simply 

. . . , cos (2?r A^5 ? - -6+ (C- 1)tt/4) , , 

c r (L,e,e) w > exp(-e) v . 3 vs 1 ' 1 30 

U VI 2 -* sin (2^.(0)1 1 ^ 

where ( = sign (2 — ft sin (27rgj(/))) and 

3/ 2 / 1 k k 2 
Sj = T + 3 '~ 2 - 3 1 + ^2 sin + cos2 ( 2 ^(0)- (31) 

In figure 2 (b) we plot the difference between the exact expression (|23|) and the semiclas- 
sical approximation (|30|) . 

For larger k, keeping terms only up to second order in the approximation (f25j) fails in 
two ways. First, 2 — k sin(27rgj(/)) vanishes for some values of / and higher terms of the 
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expansion of the action need to be taken into account. It is clear in Figure 3, where the 
predicted peak-centres (obtained by solving (|27jl and (|28j)) are plotted, that for k around 
3.0 new /-trajectories contribute to the autocorrelation with a peak around / = 0.5. 
Although these /-caustics are an additional complication, they are in fact unimportant 
if the object is to compute the contribution of bifurcating periodic orbits which takes 
place for / around zero (recall that this is the region that is semiclassically interesting for 
correlation length-scales that vanish with h). 

Second, when / — > 0, 2 — Ksin(27rgj(0)) = at periodic orbit bifurcations. Then again 
the quadratic term in (|25|) vanishes, and the fringe structure comes from the cubic term. 
It thus has a y-length-scale of the order of iV _1//3 . In the language of Section 2, this 
corresponds to a codimension-one bifurcation of a periodic orbit with r = 1 (a tangent 
bifurcation). The j = and j = 1 terms each correspond to a single unstable fixed 
point if the condition 2 — k sin(27rgj(/)) ^ is satisfied. Note that at both /-caustics and 
tangent bifucations, the semiclassical asymptotics is determined by the cubic term in the 
expansion of the action; hence both contributions are of the same order in h. 

Figure 4 provides a panorama of the situation. Here, the absolute value of the Gaussian 
smoothed autocorrelation function c is plotted. The result should be compared with the 
predicted peak-centres shown in Figure 3. The peak coming from the bifurcating periodic 
orbits is around k = 6 and / = 0. The other similar peaks are all associated with the 
/-caustics. It is clear that their amplitudes are much larger than those of the peaks 
associated with isolated /-trajectories. 

For the sine perturbation, the first bifurcation occurs when k = k* = 5.943388. At 
this parameter value, two new degenerate solutions of (|24^) appear for both j = and 
j = 1, corresponding in each case to the birth of a pair of fixed points, one stable and the 
other unstable. For / = the contributions of these bifurcating orbits when k = k* are: 

, „ , ^ N 1 ^ /4tt\ 1/3 . . cos (2-KiNSj - 9 - tt/4) , x 

^.^Ere^S?) eXP< - £) I cos (2^/3 < 32 > 

where Sj = (S(q + j) + if>(g_ -))/2 is the average of the actions of the two bifurcating fixed 
points at g +J and q -. More generally, the semiclassical contributions of the real orbits 
plus the complex orbits to the autocorrelation function give: 

* cos (27rJVSi - + (a - 1)tt/4) 

c r (L,9,e)^y exp(-e) , — + 

U v /|2- KS in(2 7 rg J (/))| 

Re ^ exp (-e)N 1 ^ L V 1 " V 2 exp (2mNSj - %9 - itt/4). (33) 



3=0 



Here a = (q+ - g_)/2 and |a] = 3| AS\/(2a 3 ), where AS = (S(q+) - %-))/2. Figure 5 
(a) shows the autocorrelation function for k = k*. In Figure 5 (b) the deviation of the 
semiclassical approximation (|33|) from the exact expression (|23|) is plotted. 
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3.2 Cosine perturbation v — 1 



Setting v = 1 in (fT8|) and the following equations allows us now to discuss a second 
order bifurcation. (Such bifurcations occur for the second iterate of the map with the 
sine perturbation, and so quantum mechanically influence U 2 (Q,Q). In the case of the 
cosine perturbation, they occur for fixed points of the map, and so influence U(Q,Q). 
This makes them easier to study numerically (Mende (1999).) 

As the order of the first non- vanishing derivative of the action S p (q) determines the 
order of the bifurcation, one can see that there are no higher bifurcations than those of 
second order for this perturbation. For j = 0, q = and k = 2, a bifurcation of order 
two occurs. For j = 1 there is no such bifurcation: in this case only tangent bifurcations 
occur, the first one at k ~ 9.208 (Mende 1999). As the tangent bifurcation was discussed 
previously, we will concentrate on the bifurcation of order two. For this, when k < 2 there 
is one real unstable fixed point (called the central fixed point), while for k > 2 there are 
two real fixed points (satellites) in addition to the central fixed point, which is now stable. 
The scarring of the wavefunctions caused by this bifurcation is shown as a dashed line 
is Figure 6, where J2 n \^n{Q)\ 2 ^e{0 — n ) — 1 is plotted for (a) k — 1 (one real unstable 
fixed point (j = 0) at q = and another real unstable fixed point (j = 1) at q — 0.5), 
(b) k = k* = 2 (where the j = fixed point bifurcates), and (c) K — 3 (three real j = 
fixed points - unstable satellites at q = —0.23 and q = 0.23 and the central fixed point at 
q — 0, which is now stable - and the unstable j = 1 fixed point at q = 0.5). A convolution 
of the data with a normalized Gaussian of width 0.02 is also shown (bold line). 

In Figure 7, | 2^q=i U(Q, Q)\ 2 = \TrU\ 2 is plotted as a function of k and N. One sees 
the contribution of the second order bifurcation of the j = fixed point around k = 2 
and that of the tangent bifurcation of the j = I fixed point around k = 9. The amplitude 
associated with the second order bifurcation is clearly semiclassically larger. Similarly, the 
amplitude of the autocorrelation function is larger at these points (compare (pSj) below 
with l|32jl). The predicted positions of the peaks in the autocorrelation function (obtained 
by solving (|2*7j) and (|28|) . when u — 1) are shown in Figure 8. Here again, for k = the 
peaks are centred at I — = 1, I — 1/3 and I = 2/3, as previously described. As k 
increases, these peaks deviate from their original positions, but the autocorrelation result 
for isolated points, the analogue of (pT0|> . is still valid. For k around 2, where the second 
order bifurcation occurs, the approximation (|25|) breaks down and an expansion of the 
phase of (fT8|) around qj up to fourth order terms is necessary. The contribution of the 
bifurcating fixed point (j = 0) when k = k* = 2 and L = is 



where S is the action of central orbit at q = 0. More generally, for k < 9 the semiclassical 



AM/4 




1/4 



Cr(O,0,e) 



2^(3/4) 



exp (-e) cos (2mNS - - vr/4) 



(34) 
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approximation to the autocorrelation function is 



cos (2nN Sl -e + (C-1W4) 
\ ^12-^008(2^(0)1 

+ Re -N 1/2 q s (l) [j_ 1/4 (z) exp (ieri) + sign(2 - «)Ji/ 4 (^) exp (icr 2 )] j exp (-e). (35) 

Here q s (l) is the distance of the satellite points from the origin, z = nNaq^/A, where 
a = 4(-g s 2 + k(1 - cos (2irq s )))/q*, a x = -9 - tt/8 + 2nNS - irNaq*/i, and a 2 = 
-9 + 5tt/8 + 2nNS - uNaq^/A. 

In Figure 9 (a) the real part of the autocorrelation function for k = k* = 2 is plotted 
and in (b) the deviation of the semiclassical approximation (|3~5|) from the exact expression 
(|23|) is shown. In Figure 10 the absolute value of the convolution of the autocorrelation 
function (|22|) with a normalized Gaussian is plotted as a function of n and /, showing the 
consistency of the theory: higher amplitudes near /-caustics and closed agreement with 
the predicted positions of the peak-centres (Figure 8). Finally, in Figure 11 the absolute 
value of the real part of the autocorrelation function (|2~3l) at L = (convolved with a 
Gaussian in k of width 0.1) is plotted. From this one sees clearly the influence of the two 
bifurcations: the tangent bifurcation at n ~ 9.208 and the second order bifurcation at 
k = 2. 
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Figure 1: c r (L,9,e) plotted as a function of / = L/N with N = 597, k = 0.0 (i.e. no 
perturbation), e = 1.5 and 9 = ir/4 (dashed line). The bold line represents a convolution 
of the data with a normalized Gaussian of width 0.007. 
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Figure 2: (a) c r (L,9,e) plotted as a function of I = L/N for the sine perturbation with 
N = 597, k = 1.0, e = 1.5 and 6 = 7r/4 (dashed line). The bold line represents a 
convolution of the data with a normalized Gaussian of width 0.007. In (b) the difference 
between the exact expression and the semiclassical approximation (jHUl) is plotted. 
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Figure 3: Predicted peak-centres in the autocorrelation function due to /-trajectories, 
computed by solving (|27|) and |28l) . for the sine perturbation. 
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Figure 4: The absolute value of the convolution of the autocorrelation function (I22J) with 
a normalized Gaussian of width 0.01, plotted as a function of k and I = L/N for the sine 
perturbation with N = 597, e = 1.5 and 9 = ir/4. A projection is shown underneath. 
Compare the positions of the ridges with the theoretical curves shown in Figure 3. Note 
also the enhanced amplitude near the bifurcation and the /-caustics. 
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Figure 5: (a) c r (L,9,e) plotted as a function of I = L/N for the sine perturbation with 
N = 597, k = 5.943388, e = 1.5, 9 = tt/4 and 6 = vr/4 (dashed line). The bold line 
represents a convolution of the data with a normalized Gaussian of width 0.007. In (b) 
the difference between the exact expression and the semiclassical approximation (|3T>|) 
is plotted. 
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Figure 6: |^„(Q)| 2 5 e (# - n ) - 1 with £ = 2.2, = tt/4 and N = 1597 (dashed line) 
for the cosine perturbation. Also shown is a convolution of the data with a normalized 
Gaussian of width 0.007 (bold line). In (a) k = 1.0 and there is only one real unstable 
fixed point (the central fixed point) at the origin; in (b) k = 2.0, which is the parameter 
value at which the second order bifurcation takes place; in (c) k = 3.0, and there are two 
new real unstable fixed points at q ~ 0.23 and q ~ 0.77 in addition to the central fixed 
point, which is still located at the origin, but now is stable. 
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Figure 7: \TrU (k, N)\ 2 plotted as a function of k and N for the cosine perturbation. Note 
the two ridges: one due to the tangent bifurcation (at k ~ 9.208) and the other due to 
the second order bifurcation (at k = 2). A projection of the data is shown underneath. 
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Figure 8: Predicted peak-centres in the autocorrelation function due to /-trajectories, 
computed by solving (127}) and (jZBj) ) for the cosine perturbation. 
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Figure 9: (a) c r (L, 9, e) plotted as a function of I = L/N for the cosine perturbation with 
N = 597, k — 2, e = 1.5 and 6 = 7r/4 (dashed line). The bold line represents a convolution 
of the data with a normalized Gaussian of width 0.007. In (b) the difference between the 
exact expression and the semiclassical approximation (I3HJ) is plotted. 




Figure 10: The absolute value of the convolution of the autocorrelation function (l2~2~|) with 
a normalized Gaussian of width 0.01 plotted as a function of k and I = L/N for the cosine 
perturbation with N = 597, e = 1.5 and 9 — 7r/4. A projection is shown underneath. 
Compare the positions of the ridges with the theoretical curves shown in Figure 8. Note 
also the enhanced amplitude near the bifurcations (see also Figure 11) and the /-caustics. 
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Figure 11: A convolution of the absolute value of the real part of the autocorrelation 
function (|2Tfl) at L = with a normalized Gaussian in k of width 0.1, plotted as a function 
of k for the cosine perturbation with N = 597, e = 1.5 and 9 = it /A. Note the peaks 
associated with the tangent bifurcation (at n ~ 9.208) and the second order bifurcation 
at (k = 2). 
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